Mori. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 16 March 2009 (MN M£X style file v2.2) 

Photo-heating and the fate of hard photons during the 
reionisation of He n by quasars 
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ABSTRACT 

We use a combination of analytic and numerical arguments to consider the im- 
pact of quasar photo-heating during He n reionisation on the thermal evolution of 
the intergalactic medium (IGM). We demonstrate that rapid (Az < 0.1-0.2), strong 
(AT > 10 4 K) photo-heating is difficult to achieve across the entire IGM unless quasar 
spectra are significantly harder than implied by current observational constraints. Al- 
though filtering of intrinsic quasar radiation through dense regions in the IGM does 
increase the mean excess energy per He n photo-ionisation, it also weakens the ra- 
diation intensity and lowers the photo-ionisation rate, preventing rapid heating over 
time intervals shorter than the local photo-ionisation timescale. Moreover, the hard 
photons responsible for the strongest heating are more likely to deposit their energy 
inside dense clumps, which cool rapidly and are furthermore invisible to most obser- 
vational probes of the IGM temperature. The abundance of such clumps is, however, 
uncertain and model-dependent, leading to a fairly large uncertainty in the photo- 
heating rates. Nevertheless, although some of the IGM may be exposed to a hardened 
and weakened ionising background for long periods, most of the IGM must instead be 
reionised by the more abundant, softer photons and with accordingly modest heating 
rates (AT < 10 4 K), although localised patches of much higher temperature are still 
possible. The repeated ionisation of fossil quasar He in regions does not increase the 
net heating because the recombination times in these regions typically exceed the IGM 
cooling times and the average time lag between successive rounds of quasar activity. 
Detailed line-of-sight radiative transfer simulations confirm these expectations and 
predict a rich thermal structure in the IGM during He n reionisation. The resulting 
complex relationship between temperature and density may help resolve discrepancies 
between optically thin simulations of the Lycv forest and recent observations. 

Key words: radiative transfer - methods: numerical - intergalactic medium - quasars: 
absorption lines - cosmology: theory. 



1 INTRODUCTION 

Helium, the second most abundant chemical element in 
the Universe, is expected to be singly ionised at the 
same time as hydrogen due to its similar ionisation po- 
tential (24.6 eV), larger photo-ionisation cross-section and 
lower abundance relative to hydrogen. Even for ionising 
sources with relatively soft, stellar like spectra, the reion- 
isation of H i a nd He i in the IGM i s expected to be 
concurrent (e.g. iGiroux fc Shapiro! I1996T ). In contrast, the 
much larger ionisation threshold (54.4 eV) and smaller 
photo-ionisation cross-section of singly ionised helium re- 
quires that He n reionisation proceeds only when sources 



_pu 



with sufficiently hard spectra become numerous, 

tativ e population of massive , metal free stars ilOh et al 
20011 : IVenkatesan et al.l 120031); high redshift X-rays (| Oh 



200ll : fvenkatesan et al J 1200 J ; iRicotti fc Ostrikej |2004) or 



therm al emission from shock heated gas ( Miniati et al.l 
2004) may provide the requisite ionising photons. The 
stacked spectra of Lyman break galaxies at z ~ 3 
l|Shaplev et al.l 120031 ) also exhibit an He n recombination 
line strong enough to suggest tha t galaxies may contribut e 
significantly to He n reionisation (|Furlanetto fc Ohll2008bh . 
However, quasars are considered to be th e most likely 
candidates for completing He n reionisation (jMadau et al.l 
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19991; iMiralda-Escude et all |2000|; ISokasian et all |2002|: 



Bolton et al.l 120061 ; IFurlanetto fc Qhl l2008bl ; IPaschos et ail 
20071 ). although bright quasars are rare at high red- 
shift and unlikely to contribute substantially to either 
the H 1 or He 11 ionising photon budget at z > 5 
JDiikstra et al J 12004 iMeiksinl 120051: ISrbinoyskv fc Wvithd 



valued and per haps even inverted (7 < 1) temperature- 
density relation JGleser et alj|2005l; ITittlev fc Meiksinll2007l ; 



Jiang et al . 2008). Conse- 



120071 ; iBolton fc Haehneltll2007bl 
quently, unless He 11 reionisation was completed by an as 
yet unidentified population of sources with hard spectra at 
z > 5, its latter stages are expected to coincide with the 
observed peak in quasar activity around z ~ 2 — 3. 

There are several pieces of observational evidence 
which may support this picture. Observations of the 
He 11 Lya Gunn-Peterson trough and patchy absorption 
in the He 11 Lya forest at z — 3 indicate a situa- 
tion analogo us to that impl i ed by the H 1 opacity data 
at z ~ 6 Oakobsen et al.l 1994 ; Davidsen et al. 1996 ; 



Heap et all l200d; ISmette et all l2002l: IZheng et al 



Shull et al] 1200-4 JReimers et alJl20oS iFechner et al 



2004 



l2006h . 



How ever, like the H 1 measurements at higher redshift 
(e.g. iFan et al. 1 120061 : iBecker et al"]l2007t K the interpreta- 
tion of the He 11 data is hampered by its insensitivity 
to the He 11 fraction once the Lya absorption becomes 
saturated. Other indirect evidence includes an apparent 
boost to the IGM temp erature at z ~ 3.3 ijSchave et al] 
120001 ; iRicotti et alj|2000l ) , measurements of metal line ratios 
which indicate the ultraviolet (UV) background is harden- 
ing a round the same redshifts (ISongailalll998l ; IVladilo et al] 
[2003)5 an d the detection of a sharp dip in t he H 1 Lya for- 
est e f fective optical depth at z ~ 3.2 (JBernardi et al] 

120031 ; iFaucher-Giguere et al] 120081 ). The latter result may 
be associated with the observe d boost in the IGM tem- 
perature at the same redshifts IjTheuns et al] l2002al ). On 
the other hand, conflicting evidence e xists for all of 
these cases (e.g . IZa ldarriag a et al] l200ll ; iKim et al] |2002| ; 
iMcDonald et al]|200ll. 120051 ). Thus, the observational evi- 
dence for He 11 reionisation remains controversial. 

Nevertheless, from a theoretical standpoint it is still ex- 
pected He 11 reionisation will have a si gnificant imp act on 
the thermal state of the IGM (fMiralda-Escude fc Reeslll994l; 
Abel fc Haehneltlll999l; i Theuns et al.l l2002al ; IPaschos et"al] 



20071 ; IFurlanetto fc Ob]|2008al ) Firstly, energy injected by 
the hard sources which reionise He 11 will raise the IGM 
temperature. This is consistent with some of the obser- 
vational constraints derived from Lya forest line widths 
which imply a boo s t of AT ~ 1-2 x 1 4 K around z ~ 3 
|Schaveet al]|2000l ; IRicotti et aT1l2000l ). The magnitude of 
the increase will depend on the spectrum of the intrinsic 
ionising radiation, and thus also on the a mount it is subse- 
quen tly filtered by the intervening IGM (|Abel fc Haehneltl 
1999). Secondly, the reionisation of He 11 will alter the 
temperature-density relation of the IGM. In the opti- 
cally thin limit, the IGM will follow a tight polytropic 
temperature-density relation, T = To A 7-1 , expected to 
arise for normalised densities A = p/{p) < 10 when photo- 
heating by a spatially uniform UV ba ckground is balanced 
by cooling due to ad iabatic expansion (iHui fc Gnedinlll997l ; 
IValageas et al] 2002). However, the assumption of a spa- 
tially uniform UV background ionising an optically thin 
IGM will break down during He 11 reionisation. Instead, ra- 
diative transfer effects and/or the inhomogeneous distribu- 
tion of the ionising sources will produce a complex, multi- 



iFurlanetto fc Ob]|2008al ; IBolton et al]|2008h 

Despite all of these studies, very little work has 
been done to establish a physical argument for exactly 
how large the temperature boost in the IGM during 
He 11 reionisation will be. Furthermore, it is unclear over 
what timescale any global temperature boost should oc- 
cur. This depends on understanding how photo-heating 
in an optically thick IGM proceeds, and thus on the de- 
tails of radiative transfer through the inhomogeneous IGM. 
The observational evidence provides valuable hints, but the 
data are still too poorly constrained to distinguish be- 
tween different models in detail. A wide variety of val- 
ues for the temperature boost during He 11 reionisation, 
5 000 K < AT < 30 000 K, have therefore been adopted 
in the literature, with the timescales for the tempera- 
ture boost varying from instanta neous to periods span- 



ning several redshift units (e.g. Abel fc Haehneltl 1999; 



Theuns et al. 2002a; [inoue fc Kamavall2003l; IHui fc Haimanl 



20031 ; Furlanctto & Oh 2008a). The main reason for this 



is that photo-heating in an optically thick IGM during 
He 11 reionisation is almost always modelled in an ap- 
proximate way. Only a handful of studies have explic- 
itly quantified the expected IGM temperature boost dur- 
ing He 11 r eionisation using radiative transfer calculations. 
lAbel fc Haehneltl (| 19991 ) were the first to point out the im- 
portance of radiative transfer through an optically thick 
IGM. Using one dimensional radiative transfer simulations, 
they found the filtering of ionising radiation through the 
IGM could increase He 11 photo-heating rates by a fac- 
tor of 1.5-2.5 above the optical ly thin values, lea d ing to 
AT ~ 10 4 K. More recently, ITittlev fc Meiksinl (|2007l) 
demonstrated that the temperature boost depends sensi- 
tively on the typical spec tra of the sources w hich drive 
He 11 reionisation. Finally, IPaschos et al] l|2007l ) have pre- 
sented a detailed, three dimensional radiative transfer simu- 
lation of He 11 reionisation. They found that radiative trans- 
fer effects could boost the IGM temperature by around a 
factor of two. However, none of these studies explicitly dis- 
cuss over what timescale this heat input may reasonably 
occur, or how the temperature boost depends in detail on 
the density distribution of the intervening IGM. 

A clear understanding of photo-heating during 
He 11 reionisation is thus vital for interpreting the ther- 
mal history of the IGM. In this paper we critically re- 
examine the photo-heating of the IGM by quasars during 
He 11 reionisation. In particular, using a combination of an- 
alytical arguments and line-of-sight radiative transfer calcu- 
lations, we shall address what the likely temperature boost 
in the IGM during He 11 reionisation is, and how rapidly 
this temperature boost should occur. Furthermore, we con- 
sider what the fate of the hard photons responsible for heat- 
ing the IGM is likely to be: where do they end up in the 
IGM, and how important is the filtering of ionising radi- 
ation through the IGM in setting the temperature during 
He 11 reionisation? 

We begin in §2, where we briefly outline the relation- 
ship between He 11 photo-heating and the photo-ionisation 
timescale. This argument is key to establishing the amount 
of heat which may be injected into the IGM within a given 
time. In §3 we present a model for estimating the maximum 
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temperature achievable during He n reionisation in an opti- 
cally thick IGM. A more detailed treatment of the filtering 
of He n ionising photons through a clumpy medium is pre- 
sented in §4, and in §5 we consider the impact of He n photo- 
heating in recombining fossil He in regions. We then com- 
pare our analytical results with detailed line-of-sight radia- 
tive transfer simulations in §6, and use these results to dis- 
cuss the impact of He n reionisation on the IGM tempera- 
ture density relation in §7. We summarise our results, com- 
pare them to current observational constraints and conclude 
in §8. 

Throughout this paper we shall assume fi m = 0.26, 
n A = 0.74, Q b h 2 = 0.024, h = 0.72, as = 0.85, n s = 0.95 
and an IGM of primordia l composition with a he lium mass 
fraction of Y = 0.24 (e.g. lOlive fe Skillmanll20o3 ). The cos- 
mological para meters are consistent with the most recent 
WMAP results (|Dunklev et alj I2008T ) aside from a slightly 
larger normalisation for the power spectrum. Unless other- 
wise stated, all distances are expressed as comoving quanti- 
ties. 



2 ON He ii PHOTO-HEATING AND THE 
PHOTO-IONISATION TIMESCALE 

We begin by considering an ideal gas parcel of density 
p exposed to a spatially uniform ionising radiation field 
with specific intensity J„ = J-2i(j / /^Hi) _aa , where J_2i = 
Jhi/10~ 21 erg s" 1 cm~ 2 Hz~* sr" 1 and z/hi are the spe- 
cific intensity and frequency at the H i ionisation threshold, 
respectively, and a s is a power-law spectral index. The ther- 
mal evolution of the g as parcel in a cosmologically expanding 
mediu m is then (e.g. iHui fc Gnedinlll997l : iFurlanetto fc Ohl 
l2008ah 

dT_2/im H 2T dS Tdp 

^"lW [ ] + 3(l + *)tf + M# ' () 

where p — [(1 — Y)(l + n c /nu)+Y/4]~ 1 is the mean molecu- 
lar weight of the gas, Jif = £\ niej is the total photo-heating 
rate per unit volume summed over the species t=(Hi, He i, 
He n] with number density rii, A is the cooling rate per unit 
volume, 5 = A — 1 is the density contrast, H is the Hubble 
parameter and all other symbols have their usual meaning. 
The first term in eq. {TJ is the dominant heating term; 
during reionisation J£° 3> A at the low d ensities and tem- 
peratures typica l of most of the IGM (|Efstathioul 1 1993 : 
iKatz et al.lll996T ). The photo-heating rate for species i =[H i, 
He i, He n] is 



4-k.L 



a\,h p {y — v\)e T " du, 



(2) 



He n reionisation. The final term corresponds to the adia- 
batic cooling driven by the expansion of the Universe, and 
is the dominant cooling term at the moderate to low den- 
sities (A < 10) still expanding with the Hubble flow. The 
relevant cooling timescale for the low density, photo-ionised 
gas in the IGM is therefore always the Hubble time, £h 
JMiralda-Escude fc Reesll 19941 ). 

Let us now assume the H i and He i in the gas parcel 
has already been ionised at some earlier time, and that the 
reionisation of He n is now completed by the radiation field 
J v . We shall further assume that the gas parcel is exposed to 
this radiation field for a time t s <C £h- The temperature evo- 
lution of the gas parcel at z ~ 3 may then be approximated 
by 



~dt 



2/iTTlH 
3fcB/0 



riHcIieHcII. 



(3) 



Ignoring the collisional ionisation of He n, which is a reason- 
able approximation in the presence of large photo-ionisation 
rates, the rate of change of riHeii is 



dni 



dt 



n-eriHelliaHelll — "-HcIlfTHell + W e aHoIl) 



(4) 



where a\ are the recombination rates for species i and in 
general the photo-ionisation rate for species i —\H. i, He i, 
He n] is 



T; = 



47rJ„ 



dv. 



(•») 



Finally, assuming the case B He in recombination timescale 
t ICC ~ 1.4 x 10 9 yrs A^T/IO 4 K)°' 7 [(l + z)/4]- 3 > t s , 
eq. (T3J) reduces to nHoii — r^^jjIdriHeii/d^l, and upon sub- 
stitution in eq. (|3]) yields 



dT ^ 2fim n 
dt 31cbP 



dni 



dt 



(-E)hcII, 



(6) 



where (-E)hoII = eHoii/rHeii is the mean excess energy per 
He n photo-ionisation. In the optically thin limit (t v < 1) 
and assuming the photo-ionisation cross-sections are pro- 
portional to v~ 3 , this re duces to (-E)hoII — /i P ^Heii/(a s + 2) 
l|Abel fc Haehneltl fl999h . Eq. © thus elucidates the key 
parameters which influence the temperature of an opti- 
cally thin gas parcel during He n reionisation; (E)ndi — 
eHeii/rHoii, which depends only on the spectral shape and 
not the specific intensity of the ionising radiation, and the 
rate of change of nHoii, which does implicitly depend on the 
intensitjQ Expressed differently, \dnucn/dt\ ~ nHeii/tion, 
where tion — Thoii i s ^^e He n photo-ionisation timescale. 
If tion 3> i s , no significant heating will occur in the gas 
parcel irrespective of the average excess energy per photo- 
ionisationlj 



where t„ = ^ . rl is the optical depth to the ionising pho- 
tons summed over the species j —[R i, He i, He n] and v\ 
is the frequency at the relevant ionisation threshold. Note 
here that J v is an intrinsic spectrum and does not include 
any reprocessing by the IGM. The second term in eq. |T} de- 
scribes the adiabatic heating or cooling due to structure for- 
mation, and is small for the modest density contrasts asso- 
ciated with most of the IGM. The third term, which results 
from changes in the mean molecular weight of the gas par- 
cel, is also small compared to the photo-heating rate during 



1 In the optically thick case the ionising radiation will be hard- 
ened and weakened by the intervening gas, but the general result 
remains the same. We assume the gas parcel is optically thin in 
this instance for simplicity. 

2 On longer timescales (2> tion) ionisation equilibrium is achieved 
and recombinations become import ant. Then in general for a 
species i, dT/dt oc pai{-E)Hcll (e.g. lTheuns!l200a) . The heating 
rate is now proportional to the IGM density and is independent 
of the intensity of the ionising radiation if the ionised fraction is 
close to unity. 
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Although this argument is applied here to a single op- 
tically thin gas parcel for simplicity, it also extends more 
generally to the IGM as a whole. A substantial temperature 
boost in the IGM can only occur globally within a timescale 
t B if the He n ionising background is sufficiently hard and 
intense. If the condition ti on < t s is not met everywhere, the 
He ii ionisation rate is too low to significantly photo-heat 
the entire IGM. Instead, temperature boosts are possible 
only in the vicinity of quasars where the local He n photo- 
ionisation rate is large. Alternatively, a large, global IGM 
temperature boost over longer timescales is still possible if 
the He n ionising background is weak, yet hard and persis- 
tent. Crucially, these scenarios depend on exactly how many 
hard photons are produced by quasars over time and the 
impact of the subsequent transfer of radiation though the 
IGM on the intrinsic quasar emission. Bearing these points 
in mind, we now proceed to discuss the expected tempera- 
ture of the IGM during He n reionisation in more detail. 



3 THE TEMPERATURE IN AN OPTICALLY 
THICK IGM DURING He n REIONISATION 

3.1 Radiative transfer and the IGM temperature 

The temperature boost in the IGM during He n reionisation 
is approximately (JFurlanetto fc Ohll2008al ) 



AT 



: 0.035/hcii (-^) (E) Hell , 



(7) 



where /hgII = n.Heii/fiHe is the He n fraction in the IGM 
when He n photo-heating commences. This assumes the 
excess energy per He n photo-ionisation is shared equally 
among all the baryons in the IGM through Coulomb inter- 
actions. In the optically thin limit, (E)jiell — ft P ^Heii/(a s + 
2) = 15.5 eV for a s = 1.5, giving AT ~ 4200 K. 

Radiative transfer effects during He n reionisation can 
substantially increase this temperature boost. An optically 
thick IGM wil l alter the spectral shap e of any incident ionis- 
ing radiation ( Haard t fc Madaull 19961 ) . In a uniform medium 
the ionising photon mean free path A„ ex i/ 3 ; high energy ion- 
ising photons have a strongly suppressed photo-ionisation 
cross-section, where a v oc v~ 3 . In an optically thick IGM 
the average excess energy per He n photo-ionisation is there- 
fore pro gressively larger the fur t her from the He n i onising 
source (JAbel fc Haehnerd 1 19991 ; iBolton et all I2004T I. This 
hardened, filtered ionising radiation can then impart a sig- 
nificant temperature boost to the IGM compared to the op- 
tically thin case. In the optically thick limit, where every 
He ii ionising photon is eventually absorbed and one may 
ignore the v~ 3 weighting from the photo-ionisation cross- 
section, (-E)Heii — h p UHeii/(a B — 1) = 108.8 eV for a B — 1.5, 
leading to AT ~ 30 000 K. Thus, a conventional interpreta- 
tion is that any large, global boost to the IGM temperature 
at z ~ 3 (inferred either directly or indirectly from observa- 
tions) may be due to photo-heating rates which are amplified 
by radiative transfer effects during He n reionisation (e.g . 
Schave et ai]|200d; Imcotti et a l. 2000; Thcu ns et al"1l2002al ; 
Bernardi et al. 2003; Fauchcr-Gigucrc et al. 2008). 

Recall, however, that irrespective of how large (T)hcII 
is, no significant photo-heating will occur in a patch of the 
IGM over timescales shorter than the local photo-ionisation 
timescale, tion — ^ncii- Ionising radiation is significantly 



weakened as well as hardened at progressively larger dis- 
tances from an ionising source. This clearly has implications 
for the AT attainable over a fixed timescale. The hardest, 
most highly filtered radiation may be too weak to signif- 
icantly heat a patch of the IGM over the short timescales 
associated with quasar activity. Quantifying the effect of this 
filtering and weakening of the ionising radiation field on the 
IGM temperature evolution is non-trivial; the temperature 
will depend on the density distribution of the intervening 
IGM, the intrinsic spectra and luminosities of the ionising 
sources, their spatial distribution and their lifetimes. De- 
tailed numerical simulations of radiative transfer are there- 
fore required to fully address this issue. We may nevertheless 
develop some intuition by first using a simple toy model for 
the effect of radiative transfer on the IGM temperature. 



3.2 A simple model for estimating ionisation 
timescales and temperature boosts in an 
optically thick IGM 

Consider a quasar which emits N ionising photons per sec- 
ond above the He n ionisation threshold where 



TV 






dv. 



(8) 



and which has a b roken power law spectrum given by (e.g. 
iMadau et ai1ll999l ) 

{i/~ 03 (2500 < A < 4600 A), 
*/-°' 8 (1050 < A< 2500 A), (9) 

u- a " (A < 1050 A). 

There are a wide range of values for the extreme UV (EUV) 
spect ral index, a s , obse rved in the quasar population at z > 
0.33 (|Telfer et al.ll2002l ). We adopt a fiducial va lue of a s = 
1.5, consistent with the mean lTelfer et al.l i|2002r i obtain for 
radio-quiet quasars. However, the EUV spectral index is a 
significant uncertainty. For example. IZheng et al.l (|l997l ) find 
a somewhat softer mean EUV inde x of a s = 1.8 ma y be 
appropriate, while on the other hand lScott et al.l (120041 ) find 
Os — 0.5 at z < 1. We therefore also consider softer and 
harder quasar spectra with a s = 2.5 and a s = 0.5, broadly 
consis tent with the dispersion in a B measured bv lTelfer et al.l 
(|2002l ). 

The rate at which He n ionising photons are produced 
by such a quasar with B-band luminosity Lb (where Lb = 
vL v evaluated at 4400 A) is then 



N 



3.6 x 10 56 s" 1 a' 1 



228 A 
1050 A 



10 12 L G 



(10) 



An Lb — L* — lO 12 !/© quasar at z — 3 yields N = (3.3 x 
10 56 , 2.4xl0 55 , 3.1xl0 54 )s~ 1 fora s = (0.5, 1.5, 2.5). Note 
the quasar spectrum which is harder/softer than the fiducial 
a s — 1.5 model also emits many more/less ionising photons 
below the He n ionisation threshold for fixed Lb. Quasars 
with Lb > L* are e xpected to dom inate th e He n ionising 
emissi vity at z ~ 3 (JHopkins et al.l 1 20071 ; iFurlanetto fc Obi 
1200801 ). We adopt L B = 1O 12 L as the fiducial luminosity 
throughout this paper. 

Let us now invoke a cut-off at some arbitrary frequency 
in the quasar spectrum above the He 11 ionisation threshold, 
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Figure 1. A simple model for estimating the He II photo-ionisation timescale and hence the maximum temperature boost around a single 
quasar for an arbitrary amount of filtering of the intrinsic quasar spectrum above the He II ionisation threshold at z = 3 (see text for 
details). The solid black lines in each panel correspond to the ionisation timescale assuming a cut-off in the intrinsic quasar spectrum 
at v' = ni/n e n, where n = 1,2,..., 20, such that the IGM is optically thin to He II ionising photons with v > v' only. The horizontal 
dotted line in each panel denotes our fiducial model quasar lifetime of 10 8 yrs, while the underlying colour scale shows the logarithm of 
the expected IGM temperature increase given by eq. JT4} for the assumed cut-off frequency. The results for the fiducial quasar model, with 
Lb = lO 12 Z/0 and q s , are displayed in the left-most panel. The remaining panels, from left to right, are as for fiducial model, but with 
a s = 0.5 and a B = 2.5. 



v' > i/Heii. One may think of this as a condition which re- 
quires t„ — ► oo for v < v 1 and t v — *■ for v > v', so the IGM 
is optically thin only to photons with frequencies above the 
cut-off. This is obviously a simplification of the full radia- 
tive transfer problem; in reality the He n optical depth is a 
continuous function of frequency and the amount of atten- 
uation will depend on the density distribution of the inter- 
vening IGM. This approach is nevertheless useful because 
it crudely mimics the effect of the filtering and weakening 
of the He n ionising radiation emitted by the quasar. The 
He ii photo-ionisation rate at a comoving distance 7? from 
the quasar is 



TheII = (1 + 2) 



4tyR 2 h v v 



dv. 



(11) 



It is then straightforward to integrate this expression to ob- 
tain 



/ 228 A 



^HcII 



V 1050 A 



10 12 L 



R 



10 Mpc 



r^ 



(12) 



The photo-ionisation rate is thus strongly dependent on the 
frequency of the assumed cut-off in the ionising spectrum, 
where FhcII oc (y' /i/Heii) - " 8-3 . On the other hand, the mean 
excess energy per He 11 ionisation increases with the cut-off 
frequency, such that 



(-E)hcIi — hpV-ndi 



i^HoII 



eii V a B + 2 / 



(13) 



Assuming ii on = r HcII < t B , substituting into eq. ([7]) yields 



AT ~ 14700 K / He ii 



v /evs + 3\ 

^Hdi \a s + 2/ 



(14) 



Hence, for a quasar with a fixed luminosity and intrinsic 
spectral shape, we may estimate the maximum tempera- 
ture boost which can be achieved through photo-heating of 
He 11 if the intrinsic spectrum has been filtered by some arbi- 
trary but constant amount. Note that this simple argument 
does not include any information regarding the likelihood 
ionising radiation is filtered by the assumed amount, nor 
the timescale for which it remains filtered; to do so requires 
a detailed knowledge of the intervening IGM density distri- 
bution. Rather, this model assumes a priori that the quasar 
radiation is already filtered by the IGM as parameterised 
by the frequency cut-off u', and remains so throughout the 
lifetime of the ionising source. We will return to the issue of 
how likely ionising radiation is filtered by a certain amount 
in {H 

Finally, we must also assume a fiducial lifetime for the 
quasar. Recent numerical studies in dicate optically brigh t 
quasars have lifetimes of t B ~ 10 7 yrs (jHopkins et al.ll2005al ) ■ 
while a wide variety of observational d ata suggest q uasar 
lifetimes vary from t s ~ 10 6 - 10 8 yrs (|Martinil 120041 ). We 
adopt the upper limit of t B — 10 8 yr - which conservatively 
maximises the effect of photo-heating - as the fiducial life- 
time, although we will discuss the effect of shorter quasar 
lifetimes later. 



3.3 The maximum IGM temperature around 

individual quasars during He 11 reionisation 

This toy model is illustrated in Fig. [T] The photo-ionisation 



timescale tic 



computed from eq. (|12[) , is displayed 



as a function of R and v' for three different quasar models at 
z = 3. In all panels the black lines correspond to the ionisa- 
tion timescale assuming v' = nviieii, where n = 1, 2, . . . , 20, 
and the horizontal dotted line denotes our fiducial model 
quasar lifetime of 10 8 yrs. The underlying colour scale shows 
the logarithm of the expected temperature increase calcu- 
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lated using eq. (|14p . The unshaded regions in the lower right 
of each panel lie below the minimum ionisation timescale 
corresponding to the unfiltered quasar radiation field, where 

v' = Z-'HoII- 

The left-most panel Fig.[T]displays the photo-ionisation 
timescales for our fiducial quasar model (Lb = 10 L®, 
a B = 1.5). Recalling that we require £i on < t B for sig- 
nificant heating to occur, we see a temperature increase 
AT ~ 30 000 K (log AT ~ 4.5), expected in the optically 
thick limit where every ionising photon is absorbed so that 
(T)hoII — h p iMeii/{ct s — 1) = 108.8 eV, is only possible very 
close to the quasar if t B < 10 8 yrs (or equivalently Az < 0.13 
at z = 3). Ionising radiation which has been hardened suf- 
ficiently to produce such a large temperature boost is too 
weak to photo- ionise the He n at R > 10 Mpc. A hard, in- 
tense ionising radiation field is in any case unlikely to be 
maintained in such close proximity to the q uasar for long 
periods. As noted by I Abel fc Haehneltl (|l999l ) (and see also 
§6), the IGM will instead quickly become optically thin to 
all He n ionising photons, thereby limiting the temperature 
boost within this region to a substantially smaller value. 
Over a quasar lifetime of 10 8 yrs, the temperature boost is 
restricted to AT < 10 4 K at ~ 40 Mpc from the quasar, and 
even less at larger distances. Note again, however, that the 
exact amount of filtering will depend sensitively on the in- 
tervening density distribution of the IGM. In addition, we 
stress these temperatures represent the maximum possible 
values within a fixed volume around a quasar, since all pho- 
tons with v < u' are neglected in eq. (I14|l . Although these 
low energy photons do increase the total ionisation rate, they 
decrease the average heating. These estimates also ignore the 
effect of any cooling on the resulting IGM temperature. 

Greater temperatures boosts can be achieved out to 
larger distances around a quasar for either a harder intrin- 
sic spectrum, a brighter luminosity, or both. The former 
increases the average excess energy per photo-ionisation, 
while the latter decreases the photo-ionisation timescale. 
The central panel in Fig. [T] shows the results for the model 
with Lb = 10 12 Lq and a B — 0.5. In addition to having 
a harder spectrum, this quasar model produces around 14 
times the number of He n ionising photons in our fiducial 
model. On the other hand, a softer, fainter quasar spectrum 
with Lb = 1O 12 L0 and a B = 2.5, as shown in the right-most 
panel of Fig. [T] yields much smaller temperatures for a given 
amount of filtering. 

This simple model clearly illustrates the difficulty of at- 
taining substantial temperature boosts over large volumes 
in the IGM within short timescales. For typical quasar EUV 
spectral indices and luminosities, except for regions in close 
proximity to the quasar, there are simply not enough hard 
photons to raise the IGM temperature substantially above 
the value expected in the optically thin limit. Furthermore, 
in this example we have assumed /hcII = 1 and t B = 10 8 yrs, 
which is aroun d the upper limit for plausible quasar lifetimes 
i|Martinill2004r ) . If the lifetimes of o ptically bright quasars are 
significantly shorter than 10 8 yrs (JHopkins et al.ll2005bl ). or 
if a patch of the IGM has /hcII < 1 as might be expected 
if fossil He in region s are common prior to the completion 
of He ii reionisation (|Furlanetto et al.l 120081 ). then the tem- 
perature boost around a single quasar will be accordingly 
smaller. Thus, for our fiducial quasar model a boost to the 
IGM temperature on timescales Az<0.13 (t B < 10 8 yrs at 



z ~ 3) in excess of ~ 10 4 K appears unlikely over scales 
of R > 40 Mpc. Given th at the average separation between 
L„ quasars is ~ 100 Mpc (|Furlanetto fc Ohll2008bT ). this im- 
plies that a rapid global IGM temperature boost in excess 
of AT ~ 10 4 K is difficult to achieve. 

One important caveat, however, is that we have applied 
the above argument to the photo-heating around individual 
quasars only. Over longer timescales, the integrated effect 
of photons emitted by distant quasars with highly filtered 
spectra could produce a more gradual yet larger IGM tem- 
perature boost. However, we shall see it is quite difficult 
to reionise the IGM with only hard photons, so the actual 
temperature boost will likely be even smaller than these es- 
timates. We address this issue in more detail in S|4] 



4 THE FILTERING OF He n IONISING 
RADIATION THROUGH AN 
INHOMOGENEOUS IGM 

Thus far we have have focused on the expected temperature 
boost in the IGM without explicitly quantifying the amount 
of filtering which occurs through an optically thick, inhomo- 
geneous IGM. Furthermore, we have not considered where in 
the IGM these hard photons are subsequently absorbed. In 
this section, we consider some simple analytical arguments 
for the amount of filtering expected, and the fate of hard 
photons. While less rigorous, these serve to motivate the de- 
tailed numerical calculations that follow later, and promote 
physical intuition. 

Specifically, we shall address two key arguments for sub- 
stantial heating during He n reionisation, which we briefly 
restate here. Firstly, we know from the lifting of the 
He ii Gunn-Peterson trough at z ~ 2.8 along multiple lines- 
of-sight that most of the H e n in the Univer s e has been 
reion i sed by this time ( e.g. iHeap et ahl |2000| ; IShull et aU 
120041 ; iFechner et al.ll20o"6r ). Thus, all that matters is the av- 
erage heat input per photo-ionisation, which as discussed 
earlier is simply (L)hcII — h p i>H B u/(a B — 1) for an optically 
thick IGM. This has to be true from a simple energy conser- 
vation standpoint, and implies significant heating. Secondly, 
the spacing between quasars is much larger than the spac- 
ing between H e 11 Lyman limit systems, which dominate the 
He 11 opacity (|Abel fc Haehneltlfl999l ; iMiralda-Escude et all 
120001 ). Thus, most of the IGM is exposed to filtered, hard- 
ened radiation. While such radiation is undoubtedly weak- 
ened, even if the photo-ionisation rate is low over long peri- 
ods of time it may again still be possible to achieve substan- 
tial heating from a weak but persistent and hard ionising 
background. Thus, while we have so far demonstrated that 
it is difficult to achieve a large temperature boost globally in 
the IGM over short timescales, it still appears to be possible 
over longer periods of time. 

While the latter scenario is certainly feasible, there are 
a number of uncertainties and caveats, which have not been 
emphasised in the literature. In particular: (i) while the aver- 
age energy per ionisation in an optically thick IGM is indeed 
(L)hoII — hpVudi I (a s — 1), dense regions preferentially ab- 
sorb hard photons. The energy per photo-ionisation in the 
low density regions may therefore be systematically less than 
the average, since the average heat input in the observable 
IGM depends on the relative collective opacity of optically 
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thin and thick systems. Note that virtually all statistical 
measures of temperature in the H i Lya forest are sensitive 
to the optically thin, underdense to mildly overdense regions. 
As an example, the low column density Lya forest (where 
Doppler line widths are directly m easurable) co rresponds to 
regions of fairly low overdensity i|Schavell200ll ). (ii) There 
are many generations of quasars. While at any given time a 
patch of the IGM may be exposed only to filtered radiation, 
over many generations a nearby quasar will eventually be 
able to directly ionise that patch with only modestly filtered 
radiation. Reionising the IGM exclusively via soft photons 
does not require an unrealistic photon budget, even with a 
substantial population of Lyman limit systems. 



4.1 The fate of the hardest photons 

Why should hard photons be preferentially absorbed in 
dense systems? We can gain some physical insight from a 
simple toy model. Let us consider a uniform IGM filled 
with a single population of absorbers of fixed mass and 
some overdensity A. For this population, the radial size 
of the absorbers R <x A -1 ' 3 , and hence the optical depth 
r oc JVhcII oc A 2 ' 3 . The effective optical depth due to the 
absorbers has two interesting limits. When A is low so that 
the absorbers are optically thin, r e ff = ^2 T% = N{t), where 
N oc a- a bsort>or oc A" 2 ' 3 is the mean number of absorbers 
encountered along a line-of-sight, (r) is the mean optical 
depth of a single absorber, and absorber is the absorber 
cross-section. In this limit, the effective optical depth is in- 
dependent of the overdensity A of absorbers: (r) ex A 2//3 , 
but for increasing overdensity their cross-section falls, and 
thus fewer absorbers, N oc A~ 2 ' 3 , are encountered along 
a line-of-sight. On the other hand, once the absorbers be- 
come optically thick, their contribution to the opacity sat- 
urates, and r cff = TV = X§g& lhm /{T) oc A" 2/3 , where 
^Hcii rbcr — Tabsorbcr/ViGM is the total mass fraction of 



He ii in the absorbers. This implies that 



^absorbc 
TIGM 



X absorber 
Hoi I 



<1,<T»' 



(15) 



Thus, on average fewer photons are absorbed in optically 
thick absorbers, since only the photosphere [t v ~ 1) con- 
tributes to the opacity. For hard photons the absorbers are 
optically thin, while for soft photons the absorbers are op- 
tically thick. Thus, a greater fraction of hard photons are 
absorbed by dense absorbers, rather than the IGM, com- 
pared to soft photons. Similarly, in the interstellar medium 
of galaxies, continuum photons can have a higher proba- 
bility of being absorbed by dusty atomic/molecular clouds 
(which are optically thin at these frequencies), compared 
to Lya photons, for which the clouds are optically thick 
i|Hansen fc OhlbOOrj )! 3 ! 

From the above, it is clear that the relative absorp- 
tion of hard and soft photons depends on the abundance 



3 Another instance of this saturation of the opacity of optically 
thick absorbers comes from the well known fact that the mean 
transmission of Lya photons is always larger in a clumpy IGM 
compared to a uniform IGM, so that the transmission is always 
dominated by underdense voids (e.g. iFan et al J 12002ft . All this 
stems from the triangle inequality, (exp [— r]} < exp [— (t)]. 



of optically thick systems, with column densities which de- 
pend on frequency as TVhcII = 1/ov oc v s . In particular, the 
fate of hard photons depends strongly on the abundance of 
high column density systems, which are optically thick to 
photons just above the Lyman edge, but optically thin to 
hard photons. Unfortunately, unlike the hydrogen Lya for- 
est, it is not possible to directly observe the column den- 
sity distribution /(NhcII, z) of the He n forest in the UV. 
One is therefore forced to theoretical model s of /(Nh c ii,z) 
|Haardt fc Madaul Il99fj ; IFardal et all 1 19981 ). There are a 
number of subtle modifications which need to be made to 
such models while He n reionisation is in progress, and we 
present detailed calculations elsewhere (Oh et al. 2009, in 
prep.). 

For now, we merely remark that /(NhcII, z) can be sig- 
nificantly model dependent, and point out the implications 
of these uncertainties. The effective optical depth experi- 
enced by a photon which has frequency Vq at redshift zq be- 
tween red shifts (zq,zi) due to discrete, Poisson dis tributed 
clouds is l|Paresce et al.lll980l ; IZuo fc Phinnevll 19931 ) 

T ff (^o, 2:0,2:1) = / dz dNmf(N m ,z) [l - e~ T "] , (16) 



where t v ~ Nm(a™ + r]a^ cU ), the frequency v — 
fo(l + z)/(l + 20), and 77 = Nncii/Nm- Note it is 
generally assumed that the H e 1 opacity is negligible 
(iMiralda-Escude fc Ostrikerl 1 19921 ; lHaardt fc Madaul Il996l ; 
IFardal et al.lll998l ). although there may be regimes where 



this is not justified. While the H 1 column density distri- 
bution /(7Vhi,2) is observationally determined as a func- 
tion of redshift, 77 is not. Thus, above 4 Ry, T c ff depends 
strongly on the unknown function t;(7Vhi, z). Given a model 
for T](Nm, z) - which usually depends on Thi, THeii and tih 
- and an assumed ionising source population of quasars 
and galaxies, Fhi and FhcII can be self-consistently com- 
puted as a function of red shift. The coupled equ ations 
have to be solved iteratively IjHaardt fc Madaulll996r ). since 
the opacity depends on the ionising background and vice- 
versa. While Thi can be compared to observations of the 
Lya forest, ThcII cannot above 2 ~ 3, when absorption 
in the He 11 forest saturates. This is unfortunate: the mean 
free path of He 11 ionising photons and hence FhoII evolves 
strongly during He 11 reionisation, a fact which can cause 
considerable uncertainty in the redshift evolution of r\. 

Let us consider the qualitative behav iour of ri(Nm,z) 
JHaardt fc Madaulll99rj : IFardal et al.lll998l ). In the optically 
thin limit, 77 = {fnciinnc)/{fmnn) « (l/12)(r H i/r H oii), 
where /; is the ionisation fraction and Fi is the ionising 
background for a given ion. In this case, ry is a constant inde- 
pendent of Nhi- However, this is no longer true at higher col- 
umn densities. At some point, He 11 beconres self-shielding 
and He hi recombines, causing r\ to rise above the optically 
thin value. At even higher column densities, H 11 recombines 
and r? oc iV - / plummets drastically. This modulation of -q 
with Nhi has been modelled by radiative transfer calcu- 
lations of slabs illumina ted by ionising radiation on both 
sides (|Haardt fc Madaul Il99r3 : IFardal et"ai]|l998l ). The re- 
sults depend on assumptions about Thi and Timi, but also 
crucially the gas density hh, which then requires an un- 
certain and model-dependent correspondence to be made 
between rh and Nhi. For instance, in the widely used mod- 
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els of lHaardt fc Madaul fl996), which assume separate fixed 
densities for Lya systems above and below the self-shielding 
limit Nhi = 1.6 x 10 1 cm" , this turnove r takes place at 
Nht = 1.2 x 10 18 cm -2 . On the other hand. IMcQuinn et al.1 
(2008), following the ansatz of ISchavel (|200l|), assume that 
nu ex N HI , and instead find that the turn-over takes place 
at JVhi = 10 15 — 10 16 cm -2 , depending on the assumed value 

Of ThoII. 

This model-dependence has little impact on 
self-consistent calculations of FhcII (at least once 
He ii reionisation is complete) . This is because He n ionising 
photons are strongly weighted toward the Lyman edge, 
and the turn-over in r\ only takes place in systems which 
are already optically thick to photons at ~ 4 Ry. For such 
photons, it does not really matter whether the optical 
depth of a system is r ~ 2 or r ~ 4; the high absorption 
probability changes only slightly. For this reason, previous 
studies correctly did not regard this model-dependence as a 
liability. On the other hand, for higher energy photons — for 
which such systems are still optically thin — the exact 
column density is important, since r tx NhoII in this regime. 
This in turn translates into considerable uncertainty 
about the frequency dependence of the opacity, and hence 
the spectral shape of He n ionising background and the 
photo-heating of the IGM. 

For instance, let us consider the mean excess energy per 
He ii photo-ionisation 



(-B)hoII = 



EHcII 



f°° dv hp(u — VH B u)c v \ucr v 



\ dv e v X v a v 



^p^HcII 

(as - i + 2) ' 



(17) 



Here 



Cj 



is the He n ionising emissivity 



[photons cm - s _1 Hz~ ], A v oc J/ is the mean free path 
for He ii ionising photons and er„ = cro(fo/z y Heii) -3 , where 
(To = 1.5 x 10 -18 cm 2 is the photo-ionisation cross-section at 
the He n ionisation threshold. For a uniform optically thick 
IGM, A„ oc v 3 and (-B)hcII — h p vu c n/(a B — 1), while for the 
optically thin case, X v is constant and comparable to the 
size of the system under consideration, yielding (-B)hcII — 
hpVHeii/(as + 2). Of course, the full inhomogeneous density 
structure of the IGM must be taken into account, using 
eq. (|16p . If we ignore the effect of He in continuum recom- 
bination radiation (which can sig nificantly soften th e spec - 
trum) , the widely used model of lHaardt fc Madaul (| 19961 ) 
yields A tx z/ 1 ' 5 , and hence (E)ueii — hpVH_eii/{a B + 0.5) (F. 
Haardt, private communication). This level of heating is in- 
termediate between the optically thick and optically thin 
case. The result is easy to understand: the turnover in 77 oc- 
curs at high values of Nm, and the assumption of constant r\ 
is a reasonable first approximation. Doing so, and assuming 
that f(Nm,z) oc N^j , we obtain: 



Tcff 



[r(/3- 



»-'i(£) 



-3(0-1) 



(18) 



where A is the comoving distance travelled by the pho- 
ton, All is the mean free path at the He 11 Lyman limit 
and the term in the brackets apply when -/V m i n = [0,0-g -1 ]. 
Note here that T (with no subscript) repr esents the gamma 
functi on. For (3 = 1.5, as assumed by lHaardt fc Madaul 
(1996) in the relevant regime iVmin = <Jq , we therefore 



obtain T e g ex v , or X v ex v 1,5 . On the other hand, 
the assumption of constant r\ is not appropriate for mod- 
els where r\ turns over at much lower column densities Nhi- 
Fo r instance, in the modi fied version of the model presented 
by IMcQuinn et al.l |2008 in which nu ex N HI , we obtain 
(E) ~ h p vu_cu/(a s - 0.5), assuming Fhcii = 10" 14 s _1 . Ac- 
counting for the thermalisation between all species, this cor- 
responds to a temperature jump of AT « 15 000K: less than 
the optically thick number of AT « 30 000K, but also signif- 
icantly greater th a n the AT « 7000 K associated with the 
lHaardt fc Madaul fl_996) model for r\. These uncertainties 
ultimately have to be observationally resolved. 

To summarise: i/high -/VhcII column density systems are 
abundant, then most high energy photons will be deposited 
there. However, the abundance of such systems-particularly 
during the course of He 11 reionisation-is far from clear, lead- 
ing to considerable uncertainty in heating rates. We will con- 
sider the effect of dense clumps on photo-heating further in 
Sj6] using our radiative transfer simulations. 



4.2 Reionising He 11 with the softer photons 

The second key argument for the importance of filtered ra- 
diation is the fact that the spacing between quasars is much 
less than the spacing b etween He 11 Lyman limit systems 
|Abel fc Haehnertlll999r ) . Let us consider an updated version 
of this argument. At z ~ 3, quasa rs with Lb ~ L* ~ 10 12 Lq 
dominate the photon budget (jFurlanetto fc Oh] l2008bl ). 
These have a space density of n q ~ 10 -6 Mpc -3 , or a mean 
separation of n" 1 ' 3 ~ 100 Mpc. Let us now reasonably sup- 
pose that these quasars can only directly ionise the IGM 
out to distances where they will on average encounter a self- 
shielded Lyman limit system. Assuming the typical size of 
an absorber with overdensity A is the local Jeans length, 
then in ionisation equili brium its H e 11 column density is 
approximately given by (|Schavell200ll ') 



AT He ii~ 1.6xl0 15 cm" 2 



T 



A 3/2 



r_i4 Vio 4 k 



9/2 



,(19) 



14 s 1 . The clump will then become 



where r_i 4 = r Hc ii/10 
self-shielded once iVHcii > 1/co = 6.7 x 10 1 7 cm" 2 . This cor- 
respon ds to a characteristic overdensity (|Furlanetto fc Ohl 
l2008bh 



^(i^rm'r 



2/3 

14- 



(20) 



iMiralda-Escude et al.l (|2000l) provide a prescription for esti- 
mating the separation between clumps of overdensity A > 
Ai calibrated to numerical simulations at z — 2-4. In their 
model, the separation between clumps of density Ai is given 
by 



A, = Ao(l + ^)[l-Fv(A0n 2/3 , 



821 



(21) 



where Fv(Ai) is the fraction of volume with A < Ai, and 
XoH(z) — 60 kms -1 is a good fit to their simulations. The 
second equality uses the fact that the lMiralda-Escude et al.l 
(J2000J) density distribution rapidly approaches an isothermal 
profile at high densities, p ex r -2 , and provides an excellent 
approximation to their expression at z — 3 (see the appendix 
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of iFurlanetto fc Ohl [20051 ). Setting R = A(A<) and using 
cqs. (112[) , (|20|l and (|21[) . we thus find that a photon at the 
He ii Lyman edge can travel a comoving distance 



Rum ~ 30 Mpc 



10 12 L 



H 



1 + 2 



-2/7 



(22) 



before it is absorbed by a Lyman limit system. We have 
ignored the weak dependence (i?n m oc y 2 / 35 ) on temper- 
ature in this expression. Since i?n m ~ 30 Mpc < n -1 ' 3 ~ 
100 Mpc, it would appear that most of the IGM is inevitably 
expo sed to filtered, hardened radiation (|Abel fc Haehneltl 
1 1999 ). While such filtered radiation is undoubtedly weak 
and results in a low photo-ionisation rate (see SJ3J), over long 
periods of time, it may still seem possible to achieve sub- 
stantial heating. 

However, there are some caveats to this argument. The 
mean separation between L* quasars is not the correct 
metric in this instance. The region over which a quasar 
ionises He n is generally substantially less than n -1 ' 3 . In- 
stead, He ii must become reionised over many generations 
of quasar for mation, since t he lifetime of quasars, t s ~ 
10 6 - 10 8 yrs (|Martini|[200i ). is much less than the Hub- 
ble time. Again for a s — 1.5, a quasar can ionise a region of 
radius 



Ri 



32 Mpc /^f 



(A/h 



U/3 



Ll 



10 12 L 



1/3 



10 8 yrs 



1/3 



(23) 



where / a b s is the fraction of photons absorbed within the 
He in region. We have again taken the upper limit for the 
quasar lifetime to be conservative with regard to the max- 
imum temperature boost. Since -Ri<-Rn m < n -1 ' 3 , most 
of the photons that ionised the IGM need not have under- 
gone significant filtering, at least until /hoII < 0.1. Note 
that this requirement becomes even stricter if the quasars 
have shorter lifetimes than 10 8 yrs|j Instead, most pho- 
tons photo-ionise and heat the region directly around their 
source quasar. Significant photo-heating by strongly filtered 
radiation is thus more likely to occur towards the end of 
He ii reionisation, once most of the He n is reionised. Even 
then, as discussed previously, much of this heating could 
potentially occur in the Lyman limit systems themselves. 
Again, it appears difficult to reionise the entire IGM exclu- 
sively with hard photons and hence achieve a large, global 
temperature boost. 

One worry with this argument might be that the pho- 
ton budget for successive generations of quasars to reionise 
the IGM is prohibitive. If every quasar only ionises a fixed, 
local region of the IGM, as implied above, then are there 
enough photons available to complete reionisation? To be 
conservative we will assume Ru m — > Ri. Then, if the filling 
fraction of ionised He in regions is Qi, on average a quasar 
will only be able to freshly reionise a new volume Vi(l — Qi) 
around it(j The rest of the photons are consumed by Lyman 



4 Note the very similar scaling with luminosity of -Ri; m oc L 2 ' 7 
and Ri oc L 1 / 3 ; hereafter we shall ignore this difference. 

5 Note that this assumes quasars are Poisson distributed. Clus- 
tering can clearly increase this penalty somewhat, since it in- 
creases the probability that a quasar is born in a previously 
ionised region. However, since the correlation length, r$ = 15.2 ± 
2.7h~ 1 Mpc, for the relatively luminous quasars observed with 
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Figure 2. Reionisation hi stories for He II assu ming an ionising 
emissivity derived from the Hop kins et al.l {2007) quasar luminos- 
ity function. The volume filling factor of ionised He III regions, 
Qi(z), is shown in the case of no photon sinks (dashed curve), re- 
combinations as in eq. (125 P with C = 1, 3 (dot-dashed curves), Ly- 
man limit systems as in eq. 1241) (thick solid curve) and both pho- 
ton sink terms with = 1 (thin solid curve). The fact that Lyman 
limit systems only allow quasars to ionise a fixed volume, regard- 
less of their previous ionisation history (as represented by eq. I25C 
does not incur a significant photon cost, and He II reionisation 
can still continue to completion. 



limit systems in the existing ionised regions. By taking this 
photon sink due to He n Lyman limit systems into account 
(and ignoring all other sinks), the filling factor evolves as 

4i=1n(l-Qi), (24) 

where hi is the rate at which He n ionising photons are pro- 
duced per helium atom. Note that this differs from the cus- 
tomary relation 

Qi = Hi - RQi, (25) 

IJMadau et ai1ll999l : lMiralda-Escude et al.ll2000h . where R(t) 
is the average recombination rate per atom in ionised re- 
gions which also includes the effective clumping factor, C, of 
ionised gas. Instead, eq. (|24[) ignores the effect of recombina- 
tions in the IGM and only considers the Lyman limit systems 
as a photon sink. It mimics the effect that C must increase 
rapidly as Qi — > 1 when den se pockets of gas must begin 
to be ionised l|Furlanetto fc Qhll20051 ); quasar He ill regions 
are sufficiently large that this could potentially happen at 
all stages of the reionisation process. Note also that since 
the recombination time in underdense regions of the IGM 
is long, fossil He in cavities which were ionised by extinct 
quas ars far outnumber activ e He in regions during reionisa- 
tion i|Furlanetto et al.ll2008l . see also §5}). The ionising emis- 
sivity required to keep the IGM ionised should thus be sig- 



the Sloan Digital Sky Survey at z > 2.9 llShen et al.l |2007[) is 
comparable to R u we expect this to be an order unity effect. 
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nificantly less than that required to initially reionise it, and 
can be maintained by the active generation of quasars. 

In Fig. [2] we show the evolution of Qi(z) according to 
these two equations, where we hav e comp uted the ionising 
emissivity using the lHopkins et al.l (|2007l ) luminosity func- 
tion. The result for no photon sinks is shown by the dashed 
curve, the dot-dashed curves include recombinations as in 
eq. (|25[1 with C — 1,3, the thick solid curve includes the sink 
due to Lyman limit systems as in eq. (|24p (but not recom- 
binations), and the result including both photon sink terms 
(recombinations and Lyman limit systems) with C — 1 is 
indicated by the thin solid curve. Note that since Qi(z) is 
simply the number of ionising photons produced minus the 
photon sink terms, it can exceed unity once reionisation is 
complete, at least in the case of eq. (|25p . 

The most important point to take away is that allowing 
quasars to only ionise a fixed volume, regardless of the pre- 
vious ionisation history of their surroundings, does not incur 
a large cost in the photon budget and He n reionisation can 
still continue to completion. Most of the IGM can there- 
fore be reionised by a direct source of soft photons; there 
are few isolated patches which have only been exposed to 
hard photons (and hence which are heated significantly). 
It is also worth noting that photo-ionisation timescales are 
then likely to be short. Unlike hydrogen reionisation, the 
size of He in regions are already controlled by He n Lyman 
limit systems before overlap, since Ri ~ _R lim . Hence, the 
mean free path and thus the ionising radiation field can 
only increase mildly - at best by a factor of 2 or 3 - during 
percolation. The strength of the radiation field after over- 
lap is constrained by transmission along all line s-of-sight in 
the H e n Gunn Peterson trough. The results of IZheng et al.l 
(J200J) suggest a fit to the He n effective optical depth of 
Teff ~ l.[(l + z)/3.8f- 5 for z < 2.8, or r Hc ii ^ l(T 14 s _1 . 
This ionising background implies ionisation timescales of 
^Heii ~ 3 x 10 6 yr <C £h- This directly contradicts the pos- 
sibility that most of the IGM is pervaded by weak, filtered 
radiation; instead, most or all of the IGM must lie within 
the proximity zone of bright source(s) which gives rise to 
a strong ionising radiation field. Such sources at a slightly 
earlier epoch would cause rapid He n ionisation. 



5 PHOTO-HEATING IN FOSSIL 
He in REGIONS 

Thus far we have considered photo-heating in regions of 
the IGM where He n is reionised only once. However, 
iFurlanetto et al.l (|2008r ) have recently pointed out that fos- 
sil He in regions will be common during the epoch of 
He n reionisation; less than 50 percent of the He in created 
by dead quasars will have enough time to recombine before 
another quasar appears. The time lag between successive 
ionisat ions of a patch of the IGM is £i ag ~ 7 x 10 7 yrs at 
z = 3 (JFurlanetto et al.ll2008l ). based on the expected num- 
ber d ensity of L > 10 115 Lq quasars from iHopkins et al.l 
(2007). Could a second round of He n reionisation in a par- 
tially recombined fossil He ill region then boost the IGM 
temperature even further, leading to temperatures in excess 
of those otherwise expected from our previous arguments? 

The case-B He in recombination timescale, t ICC = 
(QHeiiiTie) -1 , can be written as 



1.4 x 10 9 yrs / T 



10 4 K 



(26) 



whereas the time to replenish a given He n fraction, f^ u , 
is ircp ~ /hSi^so The He n fraction replenished in £i ag ~ 
7 x 10 7 yrs is thus f£? u ~ 0.05 for A = 1 and T = 10 4 K. 
Note further from eq. (0 that the temperature increase 
due to He n photo- heating is proportional to /hcII- Thus, 
since the recombination time in voids filled with He m is 
so long compared to ti &g , they will only experience a small 
fraction of the heating achieved during the first round of 
He ii reionisation. In contrast, clumps with A > 20 have 
sufficient time to recombine nearly completely before un- 
dergoing a second phase of He n reionisation, potentially al- 
lowing the He n to be fully reheated for a second time. 

However, we must also consider the cooling timescale 
to establish whether any additional temperature increase 
is actually achieved following a second He n reionisation. If 
tcoo\ < ircc , even if a clump in the fossil region has had suffi- 
cient time to fully recombine such that t rec < £i ag , the IGM 
will have cooled more quickly and any subsequent photo- 
heating by a quasar can, at most, only restore the previ- 
ous temperature in the He in region|j Fig. [3] displays the 
recombination and cooling timescales as a function of tem- 
perature for gas of primordial composition at z — 3 in photo- 
ionisation equilibrium with a UV background specified by 



Jv = J-21 ( 

*!^HI 



1 (fHI < V < fHell), 
(i/HcII < V), 



(27) 



where J_2i = 0.5 and «b = 3. The amplitude is consistent 
with measurements of the H i photo-ionisatio n rate at z = 
2 — 4 based on the observed Lya forest opacity (|Bolton et al.l 
120051 ). The solid curves show the radiative cooling timescale 



£rad — T 



tdT 
\~ctt 



3pk B T 



2fimu_{A- je) 



(28) 



at three different IGM densiti es. The cooling rate is com - 
puted using the rates listed in lBolton fc HaehneltJ (|2007al ). 
with the exception of the case-B recombi nation cooling 
rates which use the functions given by iHui fc Gnedinl 
|l997J). The dashed curves display the corresponding case- 
B He in recombination timescale. The minimum in the ra- 
diative cooling timescale at T ~ 10 4 ' 9 K results from 
He ii excitation cooling, while the high temperature tail is 
due to the inverse Compton scattering of free electrons off 
cosmic microwave background photons. The horizontal dot- 
ted line shows the adiabatic cooling timescale for the low 
density gas which is still expanding with the Hubble flow 
(A<10). For temperatures T > 10 4 K at A = 1, the re- 
combination time exceeds the adiabatic cooling timescale. 
Hence, the IGM temperature in the voids following a second 
reheating will in general be smaller than the temperature 
immediately following the initial reionisation, irrespective 
of the time lag between successive ionisations. 

At higher densities, adiabatic expansion will no longer 
dominate the cooling timescale for clumps which have al- 
ready separated from the Hubble flow. In this instance, ra- 
diative cooling becomes dominant, and the cooling timescale 



6 The exception to this is if the second quasar has an intrin- 
sic spectrum which is substantially harder than the quasar that 
originally created the fossil He III region. 
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Figure 3. Recombination and cooling timescales as a func- 
tion of temperature for gas of primordial composition at z = 
3 in photo-ionisation equilibrium with the UV background 
given by eq. J2Zj- The solid curves show the radiative cool- 
ing timescale at three different densities, from top to bottom 
A = (1, 10, 100), while the dashed curves display the corre- 
sponding case-B He III recombination timescale as labelled on 
the figure. The horizontal dotted line shows the adiabatic cooling 
timescale for low density gas still expanding with the Hubble flow 
(A<10). 



falls below the recombination timescale for T > 10 4 ' 7 K. 
Furthermore, the cooling timescale becomes comparable to 
the time lag between successive ionisations at z — 3 for 
A = 100. If icooi < tiag, even high density clumps in the 
fossil He in region which can recombine efficiently are un- 
likely to receive a substantial additional temperature boost 
following a second reheating. Note further that the radiative 
cooling timescale in overdense regions of the IGM is likely to 
be even shorter than we have assumed here. Although the 
assumption of primordial composition is a reasonable ap- 
proximation for voids in th e IGM, metal s are ubiquitous in 
higher density regions (e.g. IPettinil 12004 ) and substantially 
increas e the radiative cooling rate, depending on their abun- 
danc e (|Maio et all 120071 ; ISmith et all 120081 ; Iwiersma et all 
120081 ). This would make any net temperature increase even 
more difficult to achieve. 

Finally, at higher redshifts the time lag between succes- 
sive ionisations for a given patch of the IGM incre a ses as 
the quasar number density drops; iFurlanetto et al.l (2008) 
estimate ti ag ~ (1.4,7.2) x 10 8 yrs at z — (4,5). On the 
other hand, the recombination and He 11 collisional excita- 
tion cooling timescales at fixed density decrease as (1 + z) s . 
Hence, at higher redshifts the fossil region is even more likely 
to have sufficient time to recombine and cool before another 
quasar can reheat the IGM. We conclude that partially or 
fully recombined fossil He in regions which are reheated a 
second time are therefore unlikely to exhibit temperatures 
in excess of those expected following the initial phase of 
He 11 reionisation. Note however, heating in these regions is 



still likely to contribute significantly to the complexity of the 
relationship between temperature and density in the IGM 
during He 11 reionisation (see SJ7J). 



6 RADIATIVE TRANSFER THROUGH AN 
INHOMOGENEOUS IGM 

6.1 Numerical code 

In the previous sections we have employed analytical ar- 
guments to model the photo-heating of the IGM dur- 
ing He 11 reionisation. However, detailed radiative transfer 
simulations (ideally, with a full three-dimensional, multi- 
frequency approach) are required to fully address this 
problem. Unfortunately, such simulations are computa- 
tionally challenging, requiring very high spatial resolu- 
tion within large simulation volumes in order to accu- 
rately track the topology of He 11 reionisation, as well as 
high frequency resolution to correctly follow the thermal 
state of the IG M. Existing calculations are thus either 
one d i mensional (Haardt & Madau 19 961; Abel fc Haehnelti 



1 19991 ; iBolton et all |2004 IBolton fc Haehnelti l2007al). ne 
gleet the modelling of multip le sources ( Maselli fc Ferraral 
120051 ; iTittlev fc Meiksm1l2007l ). or resort to approximation 
schemes which do not re cover the reprocessing of the ionising 
radia tion field precisely (ISokasian et al. 2002; Pa schos et al.l 
120071 ). The latter is crucial for correctly modelling the IGM 
thermal evolution. High spatial resolution is also needed to 
resolve the overdense regions in the IGM which play an im- 
portant role in filtering the ionising radiation. Furthermore, 
large, computationally expensive simulations allow only a 
limited parameter space to be explored. 

We therefore explore the results discussed previously 
using an updated version of the one-dimensional, multi- 
fr equency photon con servi ng algorithm described and tested 
in lBolton et"aH (|2004h and lBolton fc Haehnelti (1200731 ). This 
approach is based on the monochro matic photon-con serving 
algorithm originally developed by lAbel et al.l (|l999T ). Note 
that our simulations do not model He 11 reionisation us- 
ing la rge volumes (e.g. iPaschos et alj|2007l : iMcQuinn et al.l 
2008) and will therefore not correctly capture any three di- 
mensional effects such as long range photo-heating by mul- 
tiple quasars. Instead, our models employ very high spatial 
and mass resolution and thus resolve the dense regions in 
the IGM which are important for filtering the ionising radi- 
ation. This high resolution is computationally prohibitive in 
larger, three dimensional simulations. 

The frequency integration in the updated code is now 
performed over the range 13.6 eV < h p v < 1 keV, allow- 
ing the simulations to follow the radiative transfer of soft 
X-ray and UV photons. Photons with energies in excess of 
1 keV are relatively rare and do not contribute significantly 
to photo-heating in the IGM over a typical quasar lifetime. 
Secon dary ionisations are included following the prescrip- 
tion of lRicotti et al.l ||2002|) , and t he photo-ionisation cross- 
sections from IVerner et a l. (1996) are now used instead of 
the fits from Ostcrbrockj i|1989r ). Our one-dimensional ra- 
diative transfer code does not follow the diffuse emission 
from recombination radiation, and so an option to include 
the on-the-spot approxi mation using the case -B recombina- 
tion and cooling rates of iHui fc Gnedinl |l997|) has also been 
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added. In this work we shall always use the case-B values, 
although case-A would be the more appropriate choice in 
underdense, highly ionised regions of the IGM. However, in 
this work we only consider the reionisation of an IGM with 
an initial He n fraction /hcII = 1, and so case-B is the more 
appropriate choice. 



6.2 Initial conditions 

In order to simulate the propagation of ionising radiation 
through the inhomogeneous IGM we use density distribu- 
tions drawn from a high resolution hydrodynamical sim- 
ulation run usi ng the parallel Tree-SPH code GADGET-2 
i|SpringelH2005r ). The simulation volume is a periodic cube 
27.8 Mpc in length containing 2 x 400 3 gas and dark matter 
particles. Each gas particle has a mass of 2.24 x 10 6 Mq. This 
mass resolution ade quately resolves the Lya forest at z = 3 
(|Bolton et alJl200a ). A friends-of- friends halo finding algo- 
rithm was used to identify the ten most massive haloes in the 
simulation volume at z — 3. The haloes span a mass range of 
6.67 x W^Mq < M ha i < 2.42 x 10 12 M© and lie above the 
minimu m mass expected for quasars with Lb > 1 Lp at 
z = 2-3 l|Lidz et alj|2006l ; iFurlanetto fc Ohll2008ar i. Density 
distributions were then extracted around these haloes and 
spliced with other lines-of-sight drawn randomly from the 
simulation volum e at the same redshift follow ing the proce- 
dure described in lBolton fc Haehneltl (J2007a|). This resulted 
in 30 different line-of-sight IGM density distributions, all of 
which start at the location of one of the most massive haloes 
in the simulation volume. Each complete line-of-sight is 97.2 
Mpc in length and is divided into 1024 equally spaced pix- 
els in the radiative transfer calculation, giving a cell size 
of 94.9 kpc. At mean density at z = 3 for /hcII = 1, this 
corresponds to an He n optical depth per cell of 0.1 at the 
He ii ionisation edge. 

These density distributions are used to run a total of 270 
separate line-of-sight radiative transfer simulations, divided 
into nine groups of 30. In all these simulations, we shall 
assume the H i and He i in the IGM is initially in photo- 
ionisation equilibrium with the UVB given by eq. (|27[) . The 
first group of simulations is our fiducial quasar model with 
Lb = 10 12 L Q , q s = 1.5 and t B = 10 8 yrs. Five further 
groups are run by varying one of the quasar parameters in 
the model: specifically a s = 0.5, a B — 2.5, Lb = 10 Lq, 
Lb = 1O 13 L0 and t B = 10 7 yrs. In all of these mod- 
els we assume an isothermal IGM with T = 10 4 K prior 
to He ii reionisation , consistent with the measurements of 
ISchave et al.l ((2000). The final three groups again use the 
fiducial quasar model, plus the two models with q s = 0.5 
and 2.5. However, instead of adopting an isothermal IGM 
with T = 10 4 K prior to the quasar turning o n, a power-law 
temperature-density relati on is adopted fe.g. lHui fc Gnedinl 
ll997i : IValageas et alJ[20o3 ). Specifically, we set T = LoA 7 " 1 
for A < 10, and T = TolO 7 " 1 for A > 10, with T = 10 4 K 
and 7 = 1.3, again consistent with the lSchave et al.l (2000) 
data. This initial temperature distribution resembles the re- 
sults from optically thin hydrodynamical simulations of the 
IGM l|Bolton et al.ll2008r i. 



6.3 Six example lines-of-sight 

The results for six of the simulated lines-of-sight are shown 
in detail in Figs. [4] and [5] The temperature of the IGM is 
displayed in the upper left panels, with the corresponding 
H i and He n fractions in the upper centre panels. The IGM 
density distributions along the lines-of-sight are displayed in 
the upper right panels. The lower panels display, from left 
to right, the integrated He n optical depth along the lines- 
of-sight for photons with energies 54.4 eV and 200 eV, the 
He ii photo-ionisation rates and the average excess energy 
per He n photo-ionisation. The thin curves in the latter two 
panels correspond to the relevant quantities computed in the 
optically thin limit, while the filled circles correspond to the 
distance from the quasar where the He n optical depth at 
the He n ionisation threshold first exceeds unity. The dotted 
horizontal line in the lower central panel indicates the inverse 
of the adopted quasar lifetime, t a . A selection of numerical 
results from the six lines-of-sight are summarised in Table [T] 

In each panel of Fig.|4]the solid, dashed and dash-dotted 
curves correspond to the fiducial, soft and hard quasar mod- 
els, respectively. The temperature for all three models in- 
creases somewhat within ~ 5 Mpc of the quasar. The large 
gas overdensities encountered approaching the centre of the 
haloes (typically A > 100) strongly filter the ionising radia- 
tion, exposing the nearby IGM to very hard, intense ionising 
radiation for a short period of time. The analytical estimates 
presented in Fig. [l] are consistent with such temperature 
boosts occurring close to a quasar. However, once these re- 
gions have been fully ionised, this filtering disappears, so it 
is not important for the distant IGM. 

Outside of this region, the IGM temperatures for the 
fiducial and soft quasar models are fairly modest, with max- 
imum temperatures of T — 18 700 K and 14 500 K occurring 
at R — 24 Mpc and 38 Mpc, respectively. The temperature 
in the soft model corresponds to the dense clump located 
around 38 Mpc from the quasar. This is consistent with 
many of the hardest photons being deposited in He n Lyman 
limit systems, as argued in detail in SjJ] In the hard quasar 
model a dramatic boost in the IGM temperature is ap- 
parent. The maximum temperature is T = 34 100 K at 
R = 38 Mpc, again corresponding to the overdense struc- 
ture with A ~ 50 shown in the upper right panel. The void 
which appears immediately down-stream from this clump 
is irradiated by hard, filtered yet intense ionising radia- 
tion which photo-heats the gas to a temperature substan- 
tially higher than the mildly overdense regions closer to the 
quasar. Increased temperatures in the underdense regions in 
the IGM may be required to reconcile the observed flux dis- 
tribu t ion of the Lya for est with simulations (]Becker et al.l 
I2007J ; iBolton et all 120081 ) . Note also that the presence of a 
second highly overdense clump at R ~ 47 Mpc, which has 
yet to be fully ionised by the quasar, further attenuates and 
hardens the radiation. However, photo-heating by this very 
hard but weak radiation is now less effective beyond this 
clump, where FhcII < t^ 1 . 

In general the temperatures within the ionised zone all 
increase with distance from the quasar; regions further from 
the quasar are more likely t o be ionised by hard p hotons 
with long mean free paths (|Abel fc Haehneltl |l999) . The 
temperature then falls again beyond the He in ionisation 
front, where the ionisation rate declines. The average tem- 
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Figure 4. Transfer of ionising radiation emitted by a quasar with Lb = 1O 12 L0 through a density distribution drawn from a high resolution 
hydrodynamical simulation of the IGM at z = 3. The quasar is situated in a halo on the left of each panel and has an age of t s = 10 8 yrs at 
z = 3. The gas along the line-of-sight is initially in photo-ionisation equilibrium with a spatially uniform UV background with J_2i = 0.5 
and «b = 3 between the H I and He II ionisation thresholds. Above the He II ionisation threshold the background is switched off. Three 
different EUV spectral indices are adopted, the fiducial model with ct B = 1.5 (solid curves), ct B = 0.5 (dash-dotted curves) and ct B = 2.5 
(dashed curves). Upper left: The temperature structure at z = 3 along the line-of-sight from the quasar. Upper centre: The H I (thin curves) 
and He II (thick curves) fractions along the line-of-sight. Upper right: The IGM density distribution. Lower left: The He II optical depth, 
T^ en , evaluated at the He II ionisation threshold (54.4 eV) and 200 eV. Lower centre: The He II photo-ionisation rate. The thin curves are 
the rates computed in the optically thin limit, and the dotted line corresponds to the inverse of the quasar lifetime, t B . Lower right: The 
average excess energy per He II photo-ionisation, (-E)Hell = e Hell/THell- I 11 au panels the filled circles correspond to the comoving distance 
from the quasar where 'rS ^ I eV first exceeds unity. 

Table 1. Summary of the data presented for the six example lines-of-sight in Figs. [4] and [5] In all models the initial IGM temperature- 
density relation is assumed to be isothermal with T = 10 4 K and the temperature boost, AT, is computed accordingly. The fiducial 
model corresponds to a quasar with Lb = 10 12 Lq, t B = 10 8 yrs and ct B = 1.5. From left to right, the columns list the model name, 
the mean temperature boost for all gas with /hoII < 0.1, the corresponding mean excess energy per photo-ionisation computed using 
eq. (0, the distance from the quasar where t|' c ^ i v ~ 1 and the corresponding normalised density, temperature boost, photo-ionisation 
rate and mean excess energy per photo-ionisation as a fraction of the optically thin values. The estimates for the average temperature 
boosts include the biased regions within ~ 5 Mpc of the quasar. 
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Figure 5. As for Fig. [4] except the EUV spectral index is now fixed at a B = 1.5 and three different quasar luminosities are instead adopted: 
the fiducial model with Lb = lO 12 Z/0 (solid curves), Lb = 10 13 Lq (dash-dotted curves) and Lb = 10 11 Lq (dashed curves). Note the IGM 
density distribution along this line-of-sight differs from Fig. [4] 



peratures of the reionised IGM (/hoII < 0.1) are (T) = 
(20 700, 16 200, 14 600) K for a s = (0.5, 1.5, 2.5). For an 
initial temperature of T = 10 4 K, this implies (AT) ~ 
(10 700, 6 200, 4 600) K, again consistent with the sim- 
ple analytical estimates for AT in Fig. [T] From eq. 0, 
these temperature boosts imply (-E)h«ii — (40,23,17) eV, 
corresponding to (1.8, 1.5, 1.4) times (L7)h ii evaluated in 
the optically thin limit. Excluding the biased regions 
within 5 Mpc of the quasar yields similar values; (AT) = 
(10 800, 5 900, 3 900) K. Thus on average, most of the 
photo-heating over the quasar lifetime is achieved by the 
moderately filtered radiation close to the leading edge of the 
He in ionisation front. For comparison to other estimates, 
when negle cting re-emission from th e IGM the updated cal- 
culations of lHaardt fc Madaul (| 19961 ) indicate that (-E)hcII = 
24.9 eV for a s — 1.57 within a quasar dominated UV back- 
ground model (F. Haardt, private communication). This is 
in good agreement with our results. Note, however, that in- 
cluding He in recombination emission lowers the average en- 
ergy per photo-ionisation, such that (Lv)hoII = 17.5 eV. This 
suggests that significant He m recombination emission may 
lower the expected IGM temperature boost even further by 
increasing the number of soft He n ionising photons. 

In Fig. [5] the results for the fiducial quasar model along 
a different line-of-sight are compared to two further mod- 
els with Lb = 10 11 L Q (dashed curves) and Lb = 10 13 L Q 
(dot-dashed curves). The IGM temperature again increases 
within ~ 5 Mpc of the quasar in all models. However, in the 



case of the Lb = lO n L0 quasar, there is a significant tem- 
perature boost over a region R ~ 10 Mpc. At first this may 
seem counter-intuitive; from the argument presented in $3.2\ 
one might naively expect a faint ionising source to produce 
less heating due to longer photo-ionisation timescales. How- 
ever, in this instance, the faint quasar has not had enough 
time to overcome the recombinations in its host halo, and 
as a result the radiation has been persistently filtered and 
hardened ((E)hcII ~ 38 eV). This effect is relatively com- 
mon in our simulated lines-of-sight for Lb = 10 11 Lq. Note, 
however, this effect will be less prominent or even absent 
if all the He n within ~ 5 Mpc of the quasar host halo has 
been previously photo-ionised or is instead shock-heated and 
hence collisionally ionised. 

In comparison to the line-of-sight in Fig. [4] the IGM 
density distribution in Fig. [5] is also less dense on average. 
This results in slightly less spectral hardening close to the 
quasar. As a consequence, the fiducial model in Fig. [5] ex- 
hibits lower temperatures, with a maximum of 16 500 K 
at R — 22 Mpc, while the L B = 1O 13 L model has a 
maximum of T = 18 700 K at R = 67 Mpc. The aver- 
age temperatures in regions with /hoII < 0.1 are (T) = 
(16 300, 15 400, 31 100) K for L B = (10 13 , 10 12 , 10 n )L Q . 
As before, this implies that the average excess energy per 
He ii photo-ionisation is (E)ncu — (23,20,78) eV, corre- 
sponding to (1.5, 1.3, 5.0) times (-E)Heii evaluated in the op- 
tically thin limit. Again, the first two values are consistent 
with most of the photo-heating being due to only modestly 



© 0000 RAS, MNRAS 000, 000-000 



Photo-heating during He n reionisation 15 



hardened radiation, whereas the latter model results in very 
large temperature boosts close to the quasar. Excluding the 
biased regions with 5 Mpc, the average temperatures are in- 
stead (T) = (16 200, 15 100) K for L B = (10 13 , 1O 12 )L ; 
there are no regions with /hoII < 0.1 at R > 5 Mpc in the 
Lb = 10 11 Lq line-of-sight. 



We note briefly that observations of the H i and 
He ii Lya forest are another useful diagnostic of 
He ii reionisation. Large fluctuations are observed in 
the measured column density rat i o, ?? = Nueii/Nm 
l|Zheng et al.l 12004 IShull et alj|2004 iFechner et ail l2006h . 
The average value of rj also increases towards higher 
redshifts. Direct measurements of the softness parameter, 
S = Thi/Thch <a 2 An, give S = (139±gg, 196±^°, 301±^?) 
at 2 = (2.1, 2.4, 2.8) iJBolton et al.ll200rJ l. The large error 
bars are due to the very uncertain He n effective optical 
depth measured from the He n Lya forest data. Our 
simulations are not able to follow ^-fluctuations during 
He ii reionisation in detail, since this requires a detailed 
treatment of the spatial distribution of ionising sources 
as well as radiative transfer (see e.g . iBolton et al.l 1200a ; 
IPaschos et~a?1 I2007J : iFurlanettol 120081 ). However, we may 
still estimate the average softness parameter in our sim- 
ulations along the line-of-sight from a single quasar. In 
regions with /hoII < 0.1, this yields S = (52, 153, 338) for 
the lines-of-sight with a s = (0.5, 1.5, 2.5) in Fig. [4] and 



S = (131, 237, 698) for L B = (10 13 , 10 12 , 10 11 )L Q in Fig. H 
These values are broadly consistent with the observational 
constraints, aside from the model with q s = 0.5 which 
produces a spectrum which is too hard, indicating that this 
model may be somewhat extreme. 



We have only shown the details for six of our line-of- 
sight radiative transfer calculations here, although these are 
quite typical examples from our large sample. The results 
from all the simulations will in any case be summarised in 
the next section. It is nevertheless clear from these examples 
that variations in quasar luminosities, lifetimes, and intrin- 
sic spectra, as well as the intervening IGM density distribu- 
tion, will result in a great deal of variation in the thermal 
and ionisation state of the IGM on small scales. Such de- 
tailed behaviour cannot be adequately captured by analyti- 
cal or semi-analytical arguments. They do, however, confirm 
that our general expectations hold. For typical quasar lumi- 
nosities and spectra at z ~ 3, temperature boosts during 
He ii reionisation in excess of ~ 10 4 K are unlikely to oc- 
cur rapidly (Az<0.1 - 0.2) over the entire IGM. Although 
large temperatures are indeed possible, these are most likely 
to occur either around very hard sources, or in highly over- 
dense regions of the IGM. Furthermore, the hardest photons 
tend to be preferentially deposited in the densest clumps in 
the simulations, as is evidenced by the larger temperatures 
associated with these regions. Most of the low density IGM 
is reionised by the more numerous soft photons, leading to 
temperature boosts in the range of AT = 5 000 — 10 000 K 
for a s — 1.5. Note that in all instances, there is no strong 
photo- heating wherever FhcII C t s _1 . 



7 THE IGM TEMPERATURE-DENSITY 
RELATION 

There is some recent evidence to suggest that the IGM 
temperature-density relation may be subst antially more 
compl ex th an is usually a s sumed around z ~ 3. lBecker et al.l 
l|2007h and IBolton et all (|2008h have independently noted 
that the flux distribution of the Lya forest is well repro- 
duced by a model where the voids in the IGM are signifi- 
cantly hotter than expected, perhaps due to radiative trans- 
fer effects during He n reionisation which invert (7 < 1) 
the I GM temperature-density relation. Complementary to 
this, iFurlanetto fc Ohl (|2008aT ) pointed out that the inho- 
mogeneous distribution of He 11 ionising sources will also 
complicate the IGM temperature-density relation; voids are 
generally photo-ionised and heated last and are thus hot- 
ter tow ards the end stages of He 11 reionisation. iGleser et al.l 
(2005) also demonstrated that the inhomogeneous nature 
of He 11 reionisation will blur the IGM temperature-density 
relation. We cannot address the three dimensional effects 
discussed in these latter studies with our radiative transfer 
simulations. However, we can quantify the impact filtering 
by the IGM will have on the IGM temperature-density re- 
lation. 

Contour plots of the temperature-density planes for all 
nine sets of our radiative transfer simulations are displayed 
in Fig. [6] Each panel shows the results for 30 different lines- 
of-sight. In all panels, the annotated arrows indicate the 
temperature of the IGM for all gas at mean density with 
/hoII < 0.1 (upper arrows) and /hcII > 0.9 (lower arrows). 
The dotted lines represent the initial temperature-density 
relation adopted in the models. Regions with T < 10 4 K are 
associated with He 11 which has not yet been significantly 
photo-ionised and heated. 

The upper row in Fig. [6] displays, from left to right, 
the temperature-density plane for the fiducial quasar model, 
the hard quasar model (a s = 0.5) and the soft model 
(a s = 2.5). Again note the quasar models are normalised 
to have the same rest frame B-band luminosity, such that 
N{a a = 0.5)/A>(a s = 1.5) ~ 14 and N(a 3 = 2.5)/A>(a a = 
1.5) ~ 1/8. The largest temperatures are associated with the 
hard quasar model, where additionally the most He 11 has 
been photo-ionised and heated. The middle row, from left 
to right, displays the models with Lb = 10 13 Lq, Lb = 
lO n L and t B - 10 7 yrs. Note that the L B = lO n L 
temperature-density plane exhibits the highest tempera- 
tures around mean density. This is due to photo-heating 
of the He 11 close to the quasar by intense but highly fil- 
tered radiation, as discussed in the previous section. Lastly, 
the lower row displays the same quasar models shown in 
the upper row, but now the IGM initially follows a tight 
power-law temperature density relation, T = To A 7-1 for 
A < 10 and T = T 10 7_1 for A > 10, where T = 10 4 K 
and 7 = 1.3. The temperature-density relation in an opti- 
cally thin IGM is expected to slowly asympt ote to a power- 
law w ith 7 ~ 1.6 following H 1 reionisation (|Hui fc Gnedinl 
Il997l ). This changes the absolute temperature of the under- 
dense material in the IGM following He 11 reionisation, al- 
though the temperature boost at all densities remains sim- 
ilar to the isothermal case. Thus, the IGM temperature- 
density relation prior to He 11 reionisation clearly impacts 
its post-reionisation form, although again note that the heat 
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Figure 6. Contour plots of the volume weighted temperature-density plane at z = 3 in nine different models of He II reionisation by 
quasars. Each panel shows the data from 30 different lines-of-sight through the IGM. The number density of the data points increases by an 
order of magnitude within successive contour levels. The two annotated arrows in each panel indicate the temperature of the IGM at mean 
density for regions with /h c ii < 0.1 (upper arrows) and /Hell > 0.9 (lower arrows). In all panels, the IGM is initially in photo-ionisation 
equilibrium with the UVB specified by eq. H27I I. Regions with T< 10 4 K are associated with He II which has not yet been significantly 
photo-ionised and heated. Top row: From left to right, the temperature-density plane for the fiducial quasar model (Lb = 10 12 Lq, Q s = 1.5, 
is = 10 8 yrs) and models with harder (a s = 0.5) and softer (a s = 2.5) intrinsic spectra. Middle row: From left to right, the temperature- 
density plane for the quasar models with Lb = 1O 13 L0, Lb = 1O 11 L0 and t s = 10 7 yrs. The IGM temperature in all panels in the upper 
and middle rows is initially isothermal with T = 10 4 K. Bottom row: As for the top row, except the temperature of the IGM initially 
follows a power law temperature density relation: T = ToA 7- 1 for A < 10, with T = TolO" l '~ 1 for A > 10. We assume To = 10 4 K and 
7 = 1.3. The initial IGM temperature-density relations are represented by the dotted lines in all panels. 



input into the IGM remains consistent with the analytical 
arguments presented earlier. 

The average temperatures measured at four different 
densities, log A = (—1,0,1,2), in each of the models are 
summarised in Table [2] The largest temperature boosts oc- 
cur in the densest regions in the simulations and are substan- 
tially higher than those expected in the optically thin limit 
as discussed in § [4] This suggests the hardest photons are 
preferentially deposited in these regions over the lifetime of 
the quasars. However, at mean density and below, the tem- 
perature boosts are much closer to the optically thin values, 
and only modest temperature boosts are achieved. In all 
panels there is also considerable scatter in the temperature- 
density plane. This is very different from the tight, power- 
law temperature-density relation expected for an optically 



thin IGM (e.o. lHui fc Gnedin[ll997l;IValageas et alj|2002h . A 

similar result was noted by iBolton et al.l ( 2004) following 
H i reionisation, who found that the temperature-density 
relation may even be inverted (7 < 1). However, this 
was for a single line-of-sight only and thus not necessar- 
ily rep resentative of the IGM as a whole. iTittlev fc Meiksinl 
(J2007T ) also found that radiative transfer effects can alter the 
temperature-density relation of the IGM, and will further- 
more depend on the types of sources responsible for reion- 
isation. On the other hand, using three dim ensional radia- 
tive tr ansfer simulations of He 11 reionisation. IPaschos et al.l 
1)20071 ) found the IGM temperature-density relation is not 
inverted, although the IGM temperature is still increased 
by around a factor of two at mean density. However, the 
necessarily approximate frequency integration scheme they 
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Table 2. The average temperature boost in the IGM due to He II photo-heating measured from each group of simulations displayed in 
Fig. [6] The temperature boosts are computed at four different densities, log A = (—1,0,1,2), and correspond to the temperatures in 
regions with /h c ii < 0.1 only. For comparison, the final column lists the expected temperature boost in the optically thin limit computed 
using cq. 0. 

Model (AT) [K] (AT) [K] (AT) [K] (AT) [K] AT OT [K] 

(log A = -1) (log A = 0) (log A = 1) (log A = 2) 



Fiducial 


5300 


5900 


9400 


21300 


4200 


a s = 0.5 


10400 


10600 


14900 


30100 


5900 


a s = 2.5 


3500 


3700 


4600 


16100 


3300 


L B = 1O 13 L 


6600 


6800 


9500 


21400 


4200 


L B = 10 n L Q 


5300 


7500 


15700 


26700 


4200 


t s = 10 7 yrs 


5300 


6200 


9100 


14400 


4200 


Fiducial (non-isothermal) 


5900 


5900 


7200 


15600 


4200 


a B = 0.5 (non-isothermal) 


11000 


10600 


12800 


24800 


5900 


a B = 2.5 (non-isothermal) 


4000 


3600 


2000 


9400 


3300 



use in order to run such a detailed simulation may not model 
the reprocessing of the ionising radiation field precisely. 

It is nevertheless clear from Fig.[6]that although voids in 
the IGM are significantly heated during He n reionisation, 
a well-defined, inverted temperature-density relation ap- 
pears difficult to a chieve with line-of-sig ht radiative trans- 
fer alone (see also iMcQuinn et al.l 12008 ). Empirical inter- 
pretations of the Lya forest data seem to indi cate such 
an inversion, but (as noted bv lBolton et al.ll2008h the mild 
inversion apparently required in such studies may simply 
mimic a more complex relationship between temperature 
and density in the IGM. Following He n reionisation, the 
temperature scatter in each model increases considerably, 
and over the entire IGM it may be larger still. Differences 
in quasar properties and the local density distribution of 
the IGM, as well as the thermal state of the IGM immedi- 
ately prior to He n reionisation, could enhance the overall 
scatter. The simulations also suggest the extent to which 
He ii reionisation has progressed will have an important ef- 
fect on the IGM temperature-density relation. Prior to the 
completion of He n reionisation, a two-phase medium with 
hot He in regions and cooler regions of He n should develop, 
leading to large fluctuations in the IGM temperature. Note 
that a correct treatment of the timing of He n reionisation 
at di fferent densities is also required to model this prop- 
erly IjFurlanetto fc O h 2008a), and the reionisation of par- 
tially recombined fossil He in regions will complicate the 
temperature-density rel ation yet further. On the other hand, 
studies using wavelets (JTheuns et al.ll2002bl ) have failed to 
detect spatial fluctuations in the IGM temperature on small 
scales at z ~ 3. However, these studies are implicitly pred- 
icated on fairly large (around a factor of two) temperature 
boosts; more modest boosts over small scales may be hard 
to extract from the data. 



8 CONCLUSIONS AND DISCUSSION 

In this work we have performed a detailed analysis of the 
photo-heating of the IGM during He n reionisation using a 
combination of analytical arguments and numerical radia- 
tive transfer models. Our aim was to critically examine the 
physical processes involved in modelling the fate of hard, 



He ii ionising photons in the IGM and their impact on the 
thermal evolution of the IGM. The main results of this work 
are as follows. 

(i) Filtering throug h the IGM weakens as well as hard- 
ens io nising radiation IjAbel fc Haehnehj|l999l ; iBolton et al.l 
l2004f ). However, significant photo-heating of the IGM dur- 
ing He ii reionisation can only occur over timescales com- 
parable to or shorter than the local He n photo-ionisation 
timescale, t\ ou = ThoII' irrespective of the mean excess en- 
ergy per He n photo-ionisation, (-E)hoII- Consequently, this 
limits the temperature boost, AT, attainable over a given 
timescale within a fixed volume of the IGM. For an L* 
quasar at z — 3 with a typical EUV spectra index of a> s = 1.5 
and a lifetime t B = 10 8 yrs, AT will not exceed 10 4 K at 
distances R > 40 Mpc from the quasar, and will further- 
more only do so at R < 40 Mpc if the softer photons near 
the He n ionisation edge do not contribute towards photo- 
ionisation. For an average separation between L» quasars 
of ~ 100 Mpc, this implies a global temperature boost of 
AT > 10 4 K over Az = 0.1 - 0.2 is difficult to achieve unless 
quasars with much harder EUV spectra are common. 

(ii) Hard photons are preferentially absorbed in the dens- 
est regions of the IGM where the largest temperature 
boosts will thus occur. However, current constraints on the 
IGM temperature b ased on observations of the low column 
density Lya forest (ISchave et al.1 1200(1 iRicotti et all l200d ; 
iMcDonald et alj|200ll ) are m ainly sensitive to densities close 
to the cosmic mean or below (Schayc 2001J). The fate of hard 
photons and expected heating rates thus depends on the rel- 
ative abundance of high column density iVHeii absorbers: if 
most of the heat is deposited in dense systems, we have no 
way of inferring this from the Lya forest. Unfortunately, the 
abundance of such systems is rather model-dependent, and 
varies with the unknown radiation field THeii during reioni- 
sation and the correspondence between physical density nn 
and observed column densities Nm ■ This translates into con- 
siderable uncertainties in the photo-heating rates and the 
resulting IGM temperature boost. A plausible range of val- 
ues, given model uncertainties, is AT m 7, 000 — 15, 000K, 
somewhat intermediate between the optical thin and thick 
values of AT « 4, 000K and AT w 30, 000K respectively, for 
a J„ oc v~ unprocessed source spectrum. 
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(iii) Even though the mean separation between 
He ii Lyman limit systems is comparable to the ex- 
pected size of He in regions around L* quasars, most of 
the IGM could potentially still be ionised primarily by soft 
photons, at least if quasars are randomly distributed. If 
this were the case, the average amount of heating over the 
entire IGM would be relatively modest. In general, only 
unusual regions that are far from the ionising sources will 
be exposed to a highly-filtered spectrum for a long enough 
period to heat substantially. This can only happen toward 
the end of He n reionisation, because the photo-ionisation 
timescale is so long for such heavily filtered spectra. 

(iv) The reheating of partially recombined fossil 
He in regions is unlikely to substantially increase the IGM 
temperature beyond that achieved during the first round 
of He ii reionisation. In general, recombination timescales 
in these fossil regions typically exceed cooling timescales 
and the average time lag between ionisations by successive 
generations of quasars. The net heating is therefore not 
increased. 

(v) We run a detailed set of 270 line-of-sight radiative 
transfer simulations using a variety of initial conditions 
to confirm these analytic predictions. We further demon- 
strate that the filtering of ionising radiation through the 
IGM will produce a rich thermal structure. This is likely 
to be further complicated by the inhomogeneous distribu- 
tion of the ionisi ng sources and the typ ically short life- 
times of quasars IjFurlanetto fc Ob] l2008al ) . This complex, 
multi-valued temperature-density relation may help to re- 
solve discrepancies between optically thin hydrodynamical 
simul ations of the Lyq f orest and recent detailed observa- 
tions l|Bolton et alj|2008h 

How do these conclusions compare to the current obser- 
vational constraints on the IGM thermal state at z — 3? In 
the case of the magnitude of the IGM temperature boost, 
measurements of the Doppler widths of absorption lines in 
th e Lyq forest current ly offer the best constraints. The data 
of lSchave et al.l (J2000j) are consistent with a relatively sharp 
{Az ~ 0.2), large (AT > 10 4 K) temperature jump occur- 
ring around z — 3.3 at densities around the cosmic mean. 
However, the error bars on these measurements are still 
large, and a more gradual, moderate temperature boost, as 
suggested by this w ork, is also con sistent with the data. Sim- 
ilarly, the results of lRicotti et al.l (2000) arc also consistent 
with a significant and sudden temperature boost at z ~ 3, al- 
though large statistical errors again leave a more gentle ther- 
mal ev olution possible. On the other hand, McD onald et al.l 
(2001) find no evidence for a rapid boost in the IGM temper- 
ature at 2 ~ 3. Their constraints, corrected to correspond 
to measurements at mean density, at most allow only for a 
modest evolution in the IGM temperature (AT ~ 4 000 K 
between z = 2 — 4), and are also consistent with no tem- 
perature evolution at all. A detailed analysis of a large set 
of high resolution Ly« forest spectra, perhaps using alterna- 
tive temperature sensitive statistics, is necessary to obtain 
improved constraints on the thermal state of the IGM at 
these redshifts. 

Another intriguing result which may indirectly favour 
a rapid heat input scenario, somewhat at odds with the 
results presented here, is the evidence for a sharp dip 
(Az — 0.2) of around 10 per cent in the effective optical 



depth of the Lyq forest at z = 3.2 (IBernardi et alj 120031 : 
iFaucher-Giguere et aiJioblh . lTheuns et al.l (|2002aT ) used hy- 
drodynamical simulations to interpret this feature as being 
due to a sudden increase in the IGM temperature at z ~ 3.2, 
such that AT ~ 10 4 K. However, our results suggest that it 
is unlikely this transition could occur globally over such a 
short period. Alternatively, this feature may instead be due 
to a rapid change in the hydrogen photo-ionisation rate at 
z = 3.2, although again there is no current evidence to sug- 
gest this may be the case. Further investigation into this fea- 
ture using detailed hydrodynamical simulations of the IGM 
is thus desirable. 

Our results therefore suggest a general picture for 
the thermal evolution of the IGM during He n reionisation 
by quasars where: (i) the temperature boost during 
He ii reionisation is significantly uncertain, and could be 
considerably less than has been commonly assumed in the 
optically thick heating limit, (ii) any such temperature boost 
must be achieved over a timescale longer than Az > 0.1 — 0.2 
and (iii) the resulting temperature-density relation in the 
IGM will be much more complex than the tight, power-law 
temperature density relation expected in an optically thin 
IGM. Finally, we note that these results will not apply if the 
He ii in the IGM is predominantly reionised by sources other 
than quasars. If the He n in the IGM is reionised early on by 
stellar sources, high redshift X-rays, thermal emission from 
shocked gas or some other, as yet unidentified source, then 
the thermal evolution of the IGM during He n reionisation 
may be rather different. 
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